# Lesen der Gebäudedaten
# gwr <- read.csv("https://www.web.statistik.zh.ch/ogd/daten/ressourcen/KTZH_00002022_00004064.csv") |>
# dplyr::rename_all(tolower)
gwr <- readRDS("gwr.RDS") |>
dplyr::rename_all(tolower)
#
# saveRDS(gwr, file = "gwr.RDS")
# Umwandeln in ein sf-Objekt -> Geodatenobjekt
gwr_sf <- sf::st_as_sf(gwr, coords = c("e.gebaeudekoordinate", "n.gebaeudekoordinate"), crs = 2056)
# interaktives Darstellen
mapview::mapview(head(gwr_sf, 100))
gemeinden <- sf::st_read("https://maps.zh.ch/wfs/OGDZHWFS?SERVICE=WFS&VERSION=2.0.0&Request=getfeature&TYPENAME=ms:ogd-0095_arv_basis_up_gemeinden_seen_f&outputformat=geojson") |>
dplyr::rename_all(tolower)